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ABSTRACT 



Aims. We aim to improve our picture of the low chromosphere in the quiet-Sun internetwork by investigating the intensity, horizontal 
velocity, size and lifetime variations of small bright points (BPs; diameter smaller than 0.3 arc sec) observed in the Ca ii H 3968 A 
passband along with their magnetic field parameters, derived from photospheric magnetograms. 

Methods. Several high-quality time series of disc-centre, quiet-Sun observations from the Sunrise balloon-borne solar telescope, with 
spatial resolution of around 100 km on the solar surface, have been analysed to study the dynamics of BPs observed in the Ca ii H 
passband and their dependence on the photospheric vector magnetogram signal. 

Results. Parameters such as horizontal velocity, diameter, intensity and lifetime histograms of the isolated internetwork and magnetic 
Ca II H BPs were determined. Mean values were found to be 2.2 km s \ 0.2 arcsec (~ 150 km), 1.48 {Ica) ^nd 673 sec, respectively. 
Interestingly, the brightness and the horizontal velocity of BPs are anti-correlated. Large excursions (pulses) in horizontal velocity, up 
to 15 km are present in the trajectories of most BPs. These could excite kink waves travelling into the chromosphere and possibly 
the corona, which we estimate to carry an energy flux of 310 W m~^, sufficient to heat the upper layers, although only marginally. 
Conclusions. The stable observing conditions of Sunrise and our technique for identifying and tracking BPs have allowed us to 
determine reliable parameters of these features in the internetwork. Thus we find, e.g., that they are considerably longer lived than 
previously thought. The large velocities are also reliable, and may excite kink waves. Although these wave are (marginally) energetic 
enough to heat the quiet corona, we expect a large additional contribution from larger magnetic elements populating the network and 
partly also the internetwork. 

Key words. Sun: chromosphere - Sun: photosphere - Methods: observational- Techniques: imaging spectroscopy- Techniques: 
polarimetric 
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1. Introduction point- like brigthenings as high-frequency p-modes. iLites et alJ 

^.s^ (1993) found no clear correlation between the Stokes V signal 

^ The solar chromosphere is both, heavily structured and highly *e occurrence ofK^y grains, observed in theCa n H pass- 

; dynamic. One of the important contributors to the low chromo- ^and in contrast to | Sivaraman & Livmgstoij (M who had 

. ^ sphere's structuring and dynamics in high spatial resolution im- ^""'^^"f grai ns are associate d^ a weak magnetic 

>< ages are bright points (BPs). They are also an excellent diagnos- ^0 G on average. [ Cadavid et al. | (12003) found that most 

H tic since due to their apparent lack of internal fine structure, they ^-Ime bnghtenings occur 2 minutes before or after G-band dark- 

. 5^ are ideal tracers of horizontal motions of magnetic features. Such f ^"ts in rough y the same positions in an internetwork 

BPs, however, can have a rather diverse set of physical causes, ^^ea. They concluded that the timing and coupling of these chro- 

Here we discuss particularly small BPs seen in the strong chro- mospheric brightemngs and photospheric darkenings are due to a 

mospheric lines such as Ca n H. pre-existing acoustic wave pattern with four minute periodicity. 

Acoustic waves, thought to be one of the heati ng mechanisms of They also showed that the K2V grains are not associated with any 

the quiet chromosphere from simulations (e.g.,'Carlsson & Stein" detectable magnetic flux in the simultaneously observed magne- 

[T997h an d observations (e.g., Rutten & Uitenbroek 1991; tograms. The time delay of about two minutes is roughly con- 



Kneer&vonUe xkull 1993; HofmannetalJll996; Cadavid^ ^i^^^nt the travel time of sound from the photosphere to the 
2OOI iBeck et al.. 2008 : Bello GonliteTiTa l. 2010^ annear as chromosphere (|Hoelaema.et al. 2002). 



bright, short-lived points in Ca 11 H/K filtergrams that are f". addition to the short-lived acoustic Hzv/Kzy grains, other 
sufficienfly spectrally narrow to isolate the Hzv/Kzv emission bright, point-hke structures are observed in the low solar chro- 
peak. These bright points are normally referred to as }l2v/^2v 
grains. With the help of ra d iation -hydrodynamic numerical sim- 
ulations ICarlsson & SteinI (Il997h explained the origin of these 
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Table 1. Log of observations. 



Central 


FWHM 




9 June 2009 
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12 
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12 
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8 


130 


7 








3120 


12 


150 


12 


150 


12 
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250 
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80 


12 


80 


12 


65 
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70 
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70 
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75 
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75 3 


3968 


1.8 


960 


12 


960 


12 


500 
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500 
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500 
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500 


3 


500 3 


5250.2 


0.085 


250 


33 


250 


33 














250 


33 


250 33 



* All times in UT 
** Exposure Time 
*** Cadence 



mo sphere and around the temperature minimurrQ. These include 
the br ightest parts of the (1) "reversed granulation" (iRutten et alJ 
[2001 . which appears as bright arcs all over internetwork ar- 
eas, concentrated to brighter point-like structures in some re- 
gions, (2) magnetic bright, point-like features in both network 
and internetwork areas, and (3) "persistent flashers", i.e., BPs 
which appear (when intermittent magnetic patches are squeezed 
together) and disappear (when the magnetic concentrations be- 
come less dense; as im posed by the granula r flows) during their 
lifetim e (iBrandt et alJ Il992t Ide Wiin et alJ 120051: iRutten et alJ 
[2QQ8h . 

Earlier studies of BPs in the quiet Sun mostly concentrated 
on the dynamics of network BPs. These have been well- 
studied at the atmospheric layers we sample, i.e. , from the 
upper photosphere to the lo w chromosphere (e.g., iLites et aO 
119931: (Wellsteinet"aD Il998h. They are co nnected to photo- 
spheric BPs (e.g.. iRiethmtiller et al.l 1201 Ol) with concentrated 
kilo-Gauss, nearly vertic a l magnetic elenients described by flux 
tubes (e.g., ISpruitI Il976l: ISolankil Il993l: lYeUes Chaouche et alJ 
and have been often referred to as MBPs (magn etic bright 
points; e.g., see lMuUer et al.ll2000l: iMostl et al.ll2006b or GBPs 
(when identified in the G-band passband; e.g., Berber & Title 
20011: ISteineretal.1 [2001; Sanchez A l meida et al. 2001. 20 1 ; 



Schussler et al.l l2003inShelvag et al.l l2004l: iBovelet & Wiehd 
2008*). 

The chromospheric magnetic internetwork BPs are more chal- 
lenging and less studied. The few previous in vestigations of 
magnetic internetw ork Ca ii H BP dynamics (e.g.. lAl et al.ll2002l; 
Ide Wiin et al]l2005h . although of considerable value, suff'er from 
the comparatively limited spatial resolution of the employed 
data. A complete observational picture of the dynamics of these 
bright point-like features obtained at high spatial resolution, in- 
cluding the intensity, horizontal velocity, size and lifetime dis- 
tributions, is so far unavailable. Providing such a picture for the 
smallest elements currently observable is the aim of this paper. 
Here we characterize the time-dependent properties, including 
the dynamics of the magnetic BPs observed in the low chro- 
mospheric layers sampled by the Sunrise Filter Imager (SuFI) 
Ca II H filter, centred at 3968 A. We measure their magnetic 
property with the help of simultaneous full Stokes photospheric 
observations in the Imaging Magnetograph experiment (IMaX) 



^ Note that the term "bright point" has been interchangeably used in 
the literature for chromospheric phenomena such as brightenings asso- 
ciated with magnetic features, or K- grains (e.g., Steffens et al. 1996), 
and for bright, point-like features in the other layers of the solar atmo- 
sphere, e.g., X-ray bright points in the corona, or photospheric magnetic 
bright points. Here we restrict ourselves to the chromospheric features. 



Fe I 5250.2 A passband. Several criteria are considered to define 
the particular breed of Ca ii H BPs probed in this work. 
This paper is organised as follows: in Sect. 2, we describe the 
data and their preparation prior to analysis. In Sect. 3, our meth- 
ods for studying the temporal evolution of the identified Sunrise 
Ca II H BPs, as well as their magnetic correspondence measure- 
ments are explained. The results of our investigations and related 
discussions are presented in Sect. 4, followed by the conclusions 
in Sect. 5. 



2. Observations and Data Preparation 

Several simultaneous imag e sequences obtai ned by the Sunrise 
Filter Imager (SuFI; Gand orfer et al.l 1201 ll) and the Imaging 
Magnetograph experiment (IMaX; iMartinez Pillet et aP 
[2011) on board the Sunrise balloon-borne solar tele- 
scope (ISolankietal.1120101: iBarthol etal]l201ll: iBerkefeld et"aD 
I2OIII) are used for this study. 

Table [T] summarises the SuFI and IMaX observations used in 
this study. In this table, we present exposure time and observing 
cadence of each image sequence and for all wavelengths 
recorded by Sunrise. However, only the SuFI Ca 11 H (centred at 
3968 A), SuFI CN (centred at 3880 A) and IMaX Fe I 5250.2 A 
images were used for this paper. Image sequences that were 
not observed in a certain wavelength and/or data that were not 
available (e.g., not yet reduced) at the time of this investigation, 
are indicated by dashed lines. 



2.1. Can H Images 

The primary time series in this investigation consist of SuFI im- 
ages (Hirzberger et al. 2010) taken in the Ca 11 H filter (centred 
at 3968 A with FWHM^^ 1.8 A). The series were recorded on 
2009 June 9 (at and 1 UT, both with 12 sec cadence), June 
11 (at 15, 20 and 21 UT, with 4, 8 and 7 sec cadences, respec- 
tively) and June 13 (at 1 UT, with 3 sec cadence) in quiet-Sun 
regions close to disc-centre. We employ level 3 data, i.e., images 
that have been phase diversity reconstructed with averaged wave 
front errors (cf. iHirzberger et al.ll2010ll201 ih . 
To obtain an estimate of the formation height of the radiation 
passing through the Sunrise/SuFI Ca 11 H passband, we deter- 
mine d the contribution functions (CF) by means of the ID RH- 
code (IUitenbroekl200 ih . This code solves both radiative transfer 
and statistical equilibrium equations in non-LTE for a given at- 
mospheric m odel. We used a fo ur-level plus continuum Ca 11 H 
model atom (IUitenbroekll200 ll) and three diff'erent atmospheric 
models, FALC, FALF and FALP, which respectively describe 
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[arcsec] [arcsec] [arcsec] 

Fig. 1. Example of co- spatial and co-temporal images from both SuFI and IMaX instruments on board the Sunrise balloon-borne 
solar telescope, (a) SuFI Ca ii H image, (b) Average of two line-positions at -40 mA and -h40 mA from the Fe I 5250.2 A line-core. 
Plotted is the part of the full FOV overlapping with the SuFI Ca ii H image, (c) Full field of view of the continuum position of 
the IMaX Fe I 5250.2 A line. The red dashed-line in this image illustrates the co-aligned area in common with the SuFI Ca ii H 
image. The yellow box (solid-line) includes a sample Ca ii H BP studied here; its temporal evolution and parameter fluctuations are 
presented in Fig. (Hand Fig. [51 respectively. The common coordinate system adopted for the present analysis is given below and to 
the left of panels a and b, as well as above and to the right of panel c. 



an averaged quiet- Sun area, bright regions of the qu iet-Sun and 
a typical plage region (iFontenla et al.lll993L [2QQ6b . For all at- 
mospheric models, the intensity profiles (synthetic spectra) of 
the Ca II H line, computed with the RH-code from the atomic 
model, as well as the determined contribution functions to line 
depression (e.g. j Grossmann-Doerth et a l."l988') were convolved 
with the transmission profile of the Sunrise/SuFI Ca ii H fil- 
ter. Fig. [2^ shows the Ca ii H spectra before convolution with 
the transmission profile (dashed curve) for FALC (red line), 
FALF (green line) and FALP (blue line). Fig. ^ displays the 
convolved profiles. The latter model represents higher Il2v/^2R 
emission peaks around the line-core compared to the other mod- 
els due to excess hea ting related to the magnetic field, just as 
in magnetic BPs (e. g . , ISkumanich et al.|[T975l: lAyres et alll 19861: 
ISolankietal.]ll99lh . For a comparison of Ca ii K line shapes 
in diff'erent model atmospheres see Fig. 7 of iFontenla et al.l 
(|2009). The mean intensity contrast of FALP relative to FALC 
I A3; determined from the convolved spectra) is similar to 
that of small, internetwork Sunrise Ca ii H BPs under study (see 
Sect. 14.3. it . Integration of the (filter-profile convolved) CF over 
wavelength gives the contribution of each height to formation of 
the signal seen in the Sunrise Ca ii H bandpass. Plotted in Fig. [21: 
are the computed CFs for the Sunrise filter versus geometric 
height for the three model atmospheres. The CFs have multiple 
peaks and cover a wide range of heights from below the surface 
to the mid-chromosphere. The vertical do tted lines in Fig [2^ in- 
dicate the average formation heights (cf . . iLeenaarts et al]|20 1 2) . 



lying at 437, 456 and 500 km for FALC, FALF and FALP, respec- 
tively. The sudden drop of the CF's curve for FALP at ^ 1700 km 
is due to the rapid increase in the temperature at the top of the 
chromosphere. 

The Sunrise/SuFI Ca ii H filter is narrowe r than that of Hinod e 
(centred at 3968.5 A with FWHM^ 3.0 A: lTsuneta et~ani2008h . 
so it samples higher layers in the quiet- Sun compared to the 
Hinode Ca ii H passband with a n average formation height of 
247 km. fcf. lCarlsson et al.ll2007h . 

The photon noise is the most prominent noise source in inten- 
sity and hence it is assumed as the only noise contribution in 
the measurements. The noise introduced by instrumental eff'ects 
into the Sunrise observations is much smaller and not consid- 
ered further. The photon noise is calculated by measuring the 
high frequencies derived from the Fourier power spectrum of the 
reconstructed images (after apodization). A value of 0.01 {Ica) 
was measured as the averaged photon noise in our Ca ii H im- 
ages. 



2.2. Stokes Measurements 

Full Stokes vector (I,Q,U,V) data recorded in parallel to the 
Ca II H sequences by IMaX in the magnetically sensitive pho- 
tospheric Fe I 5250.2 A line at a cadence of 33 sec were also 
analysed whenever available. The Stokes vector was measured 
in four wavelength positions in the line located at -80, -40, -h40 
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Fig. 2. (a): Intensity profiles (synthetic spectra) of Ca ii H line 
for three atmospheric models (see text). The dashed curve repre- 
sents the transmission profile of the Sunrise/SuFI Ca ii H filter. 
(b): Results of the convolution of the spectra with the transmis- 
sion profile shown in panel (a), (c): Line depression contribution 
functions for the Sunrise/SuFI Ca ii H bandpass vs. height. The 
vertical lines indicate the corresponding average heights forma- 
tion. Red, green and blue lines in all panels represent the FALC, 
FALF and FALP model atmospheres. 



and -1-80 mA from line centre and at -1-227 mA in the continuum. 
The noise level wa s 3 x 10"^ 7^ after phase-d iversity reconstruc- 
tion (V5 - 6 data :']V[artmez Pi llet et al.ll20ri ). 
In order t o study the m agnetic prop erties of the investigated BP s, 
following [Solanki et al.1 (|2Q1Q|) and llMartmez Pillet et al] (l2Qlli) . 
we averaged over Stokes V at the 4 wavelength positions inside 



the line, normalised to the local quiet-Sun continuum intensity 
(Ic). To avoid cancellation eff'ects, the sign of the two red wave- 
length points was reversed prior to averaging. This wavelength 
averaged Stokes V will be called CP in the following. 
The total linear polarisation (LP) was also computed from the 
Stokes Q and U profiles (LP = VG^TTT^), averaged over the 
four wavelength positions in the line. To do this, we first average 
the Stokes Q and U profiles over the wavelength positions with 
the same sign (i.e., the two outermost and two innermost wave- 
length positions), then compute the LP for each averaged posi- 
tion, and finally calculate the mean value of the two LPs. The CP 
and LP were determined from the non-reconstructed data (with a 
factor of 2 lower spatial resolution and, more importantly, a fac- 
tor of 3 lower noise compared to the reconstructed data), which 
were flat-fielded and corrected for instrumental efl'ects. 
The noise level for LP was computed from the IlMaX Stokes 
Q and U frames of the continuum wavelength and the LP was 
corrected for the ofl'set produced due to squaring the Stokes pro- 
files. Furthermore, to increase the signal-to-noise ratio (SNR), 
an average over a 3 x 3 pixels area centred at the corresponding 
location of the BP on the aligned IMslX images was performed 
on stokes Q, U and V profiles, at each wavelength position, prior 
to computing the CP and LP signals. The 3cr level of the Stokes 
profiles decreases due to the averaging over four wavelength po- 
sitions and it is additionally reduced due to the spatial averaging 
of the polarisation signals. The reduction in the noise by a factor 
of 3 due to the spatial smoothing over 3x3 pixels was not fully 
achieved due to coherent noise, probably caused by image jitter. 
The final computed value of the 3cr noise levels for CP and LP 
normalised to Ic is on average 1 x 10"^. 

Areas where the CP signal is greater than the computed 3cr noise 
level are defined as magnetic regions (patches). 



2.3. Image Alignment 

Since the Ca ii H images form the primary data set considered in 
this paper, the IlMaX 50'' x 50'' images (936 x 936 pixels in size 
and with an image scale of 0.0523 arcsec/ pixel) were cropped 
to the field of view (FOV) of the SuFI Ca ii H images covering 
712 X 1972 pixels (with an image scale of 0.0207 arcsec/ pixel) 
after interpolating the images to an approximately common scale 
and orientation. 

Furthermore, the images were properly aligned by applying a 
cross-correlation technique on common sets of manually se- 
lected features in two image series having the closest forma- 
tion heights: Ca ii H and the average of the -40 and -h40 mA 
line-positions of the IlMaX Stokes /. The latter corresponds to 
the inner flanks of Fe I 5250.2 A and is a proxy for the line- 
core intensity. The differences in solar position of the manu- 
ally selected features were iteratively minimized using the Powel 
minimization function. The resulting cross-aligned images show 
good alignment within the spatial resolution of the observations. 
After co-alignment, the UMaX images were given the same co- 
ordinate system as the SuFI Ca ii H images for each image se- 
quence, so that any spatial structure, e.g., a BP, is addressable by 
the same coordinates in all images. Note that the SuFI and IlMaX 
data have different cadences, hence, we selected those Ca ii H 
images whose observing times are closest to the IlMaX line-core 
observations. 

Fig. [T] shows example frames of SuFI Ca ii H (left), IlMaX line- 
core (as described above; centre) and the IMaX Stokes / contin- 
uum (right) after the co-alignment. 
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3. Analysis 

In Sect. [T] we reviewed point- like brightenings observed around 
the temperature minimum, i.e., Il2v/^2v grains, product of re- 
versed granulation, magnetic BPs and persistent flashers. They 
are in many respects similar to the BPs considered here, but 
many of the BPs in the literature also difl'er in important ways. 
In the present work, we investigate the Sunrise/SuFI Ca ii H BPs 
which (1) are located in internetwork areas, (2) are intrinsically 
magnetic (as far as we can determine), (3) are not oscillations or 
wave-like features, (4) are not the product of reversed granula- 
tion, (5) are actual bright point-like features, i.e., have a roughly 
circular shape, (6) are brighter than the mean intensity of the im- 
age and (7) do not interact with other BPs (i.e., merge or split) in 
their observed lifetimes. These conditions restrict us to internet- 
work, magnetic, non-interacting and small BPs seen in the height 
sampled by the Sunrise/SuFI Ca ii H passband. These BPs turn 
out to be the smallest Ca ii H BPs observed so far. The small 
size is partly due to the upper limit of 0.3 arc sec that we set. 
This facilitates precise tracking and avoids any ambiguities due 
to a changing shape and/or varying distribution of the brightness 
over the magnetic internetwork BPs. In order to isolate BPs with 
the above properties, we set the following criteria: 



(1) All the BPs are investigated in quiet internetwork areas 
where no significant large bright features, i.e., no obvious 
magnetic structure representing the network, are present in 
the immediate vicinity, i.e., within 3 arc sec. 

(2) The magnetic properties of the BPs are determined from the 
simultaneous photospheric magnetograms (Stokes V signals) 
obtained from the Sunrise IMaX instrument. Only BPs with 
CP above 3cr during at least half of the BPs' lifetime are 
considered. This criterion is imposed whenever simultaneous 
IMaX observations are available (see Table [B. 

(3) Criterion (2) ensures that the analysed BPs are not due 
to non-magnetic phenomena such as acoustic waves. Since 
the magnetic information is not available for the image se- 
quences observed on 11 June 2009, we per form a furthe r 
test for these time series proposed by Cadavid et al. (2003). 
They found a delay between photospheric intergranular dark- 
enings and H/K-linc brightenings. We investigate such an 
occurrence for our candidate BPs by comparing the aligned 
simultaneous SuFI images in the CN 3880 A and Ca ii H 
passbands. Thus, a BP is included only if its brightening is 
not preceded by a CN darkening. 

(4) The bright arc- shaped features caused by the reversed gran- 
ulation may in some cases produce point-like brigthenings 
and hence may lead to false BP detection. These arcs and as- 
sociated BPs are, however, short-lived. Therefore, in order to 
avoid misidentifi cation, we exclude BPs with lifetimes less 
than 80 seconds (Ide Wiin et al.ll2005h . 

(5) We set an upper limit of 0.3 arcsec, i.e 2 times the spatial 
resolution of the Sunrise data, for the BP diameter. This size 
threshold limits the identified bright features to actual point- 
like features with roughly circular shapes. 

(6) A BP is considered to be detected only if its intensity (i.e., 
contrast) in a Ca ii H image, at any given time, rises above 
the mean intensity of the entire frame ({Ica})- 

(7) Finally, BPs which merge with other features or split into 
several BPs (or magnetic elements) during their lifetimes are 
excluded. In order to sort out BPs aff'ected by such interac- 
tions an additional careful manual inspection is made of the 
temporal evolution of each BP. 



3.1. Tracking Algorithm 

Identification and tracking of Ca ii H BPs, i.e., determining 
their positions and other parameters in cons ecutive frames, is 
not straight forward as pointed out by de Wijn et all dlOOS). We 
employ a tracking algorithm b ased o n a particle tracking ap- 
proach developed by Cro cker & Gried (IT996) . This method has 
been widely used and has been confirmed as a precise approach 
in colloidal studies for identifying and tracking small roughly 
point-like features on a noisy and variable back ground (e.g., 
Croc ker & Griejll996l: IWeeks et al.ll2000l: Uenkinsa & Egeltiaal 
2008\ conditions very similar to the Ca ii H images. 
After the BPs were automatically identified in all frames, we 
again visually checked each BP in each image in order to catch 
and discard misidentifications, i.e., the interacting BPs. Then, 
starting from the frame in which a particular BP was first iden- 
tified we tracked the selected BP automatically backward and 
forward in time for as long as all criteria listed above were met. 
The brightness, size and shape of a BP change as it migrates. 
The tracking algorithm must deal with changing properties and 
distinguish the BP from noise and other transient bright, point- 
like features (e.g., caused by a passing wave) in the immediate 
surroundings of the candidate BP. 

A region slightly larger than the area that the selected BP is mov- 
ing in (i.e., a typical area of 3'' x 3'') is cropped from the whole 
field of view. This smaller subframe generally includes fewer 
non-interesting bright features that may cause difficulties to the 
tracking procedure. 

The algorithm is then initialized by determining all local in- 
tensity minima and maxima in the first subframe (the reference 
frame), followed by two steps applied to all images of the whole 
time series: (1) image restoration and (2) locating the BP and 
linking the locations into trajectories. 



3.1.1. Image Processing 

In order to facilitate precise tracking of the identified BPs, the 
image is restored to correct for inhomogeneities such as non- 
uniform background solar intensity, noise and geometric distor- 
tion due to rectangular pixels. The restoration includes a cor- 
rection for sudden brightenings in the immediate surroundings 
of the BP, e.g., caused by a passing wave. This eff'ect includes 
local intensity peaks of varying levels which in principle can- 
not be distinguished from the BP even by applying an intensity 
threshold. In this study, we follow a sophisticated method based 
on the algorithm developed by Crocker & Grier (1996) to track 
suspended particles. The algorithm restores the image by apply- 
ing a convolution kernel in which (1) the low-frequency mod- 
ulation of the background intensity with a non-uniform bright- 
ness is subtracted from the entire image after boxcar averaging 
over a circular region with a diameter of 2w -h 1 pixel, where 
2w is an even valued integer slightly larger than a single BP's 
apparent diameter in pixels, and (2) purely random noise with a 
correlation len gth of An ^ I (i.e., assum ing single pixel instru- 
mental noise; Jenkinsa & Egelhaaf 2008) is suppressed without 
blurring the image. The latter step is performed by convolving 
each image with a Gaussian kernel. In order to maintain consis- 
tency with our BP definition, we set the former size threshold 
2w slightly larger than the maximum diameter of the defined 
BPs (i.e., 0.3 arcsec) which equals 20 pixels. Assuming a fixed 
size of 2w for the BPs, we compute both these steps with the 
convolution kernel 



1 / f+f 



B 



1 



(2w + ly 



(1) 
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where and B are normalisation constants. Thus, the filtered 
image (A /) after the restoration is given by 

w 

A fix, y)= ^ A(x - /, J - j)K{i, j) , (2) 

where A(x, y) is the original image, (x, y) and (/, j) are pixel co- 
ordinates in the image and kernel, respectively. As an example, 
panel (b) in Fig. [3] shows the result of such an ima ge restoration 
applie d to the original image shown in panel (a). iBerger et^ 
([l995h restored G-band images by applying a somewhat similar 
kernel from a "blob finding" algorithm used to identify GBPs, 
followed by a further intensity enhancement of dilation and ero- 
sion processing. 

To summarise, we applied a real-space, spatial bandpass filter 
by (wavelet) convolution of the Ca ii H images with (1) a "large 
kernel" to cut off' low spatial frequencies (i.e., eliminating the 
larger scale residual noise of the bright, extended background 
structures; the second term in Eq. ([T])), and (2) a "small kernel" 
to truncate high frequency noise (the first term in Eq. ([T])); while 
retaining information of a characteristic size (2w -h 1). As a re- 
sult, we enhanced the BP's intensity profiles at the expense of 
other image components, e.g., passing waves and reversed gran- 
ulation. However, the enhanced images were only used to locate 
the BPs and to measure their sizes. They were not used in subse- 
quent measurements of the BP intensity (i.e., contrast), since the 
restoration procedure may introduce spurious high spatial fre- 
quency artifacts in intensity. The contrast values of the BPs are 
instead measured from the original Ca ii H images, by referenc- 
ing their accurate positions and sizes obtained from the restored 
images. 



0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 
[arcsec] [arcsec] 

Fig. 3. A cropped subframe including a visually selected BP 
(panel a) is restored after its background is subtracted off' and 
noise is reduced (panel b). 



3.1 .2. BP'S Location and its Trajectory 

The position of the maximum intensity in the identified (se- 
lected) BP in the reference frame is assigned the initial coor- 
dinates of the BP. These initial coordinates are then assumed to 
be at the centre of a circular window with a fixed diameter of 
20 pixels ^ 0.4 arcsec (i.e., a window slightly larger than the 
maximum allowed size of the BPs with diameter of 0.3 arcsec). 
If more than one pixel has the same maximum intensity value, 
the pixel halfway between them is selected instead. Then, (1) 
the centre of gravity of the intensity is computed for pixels in- 
side the circular window as new coordinates of the BP, (2) the 
circular centroid is moved to the new coordinates and (3) the co- 
ordinates are iteratively checked and refined if the new centre is 
found to have a larger off'set than 0.5 pixel from the current cir- 
cular window centroid. 



In the next step, the selected BP in the reference frame is com- 
pared with the automatically identified BP in the vicinity of its 
position in both the preceding and the following frames. Finally, 
the locations of a BP in all consecutive frames are linked to- 
gether to describe its trajectory. 

The displacement of a BP between consecutive already analysed 
frames can be used as an input in the tracking algorithm and is 
useful when a BP disappears and re-appears as happens for the 
so-called "persistent ffashers". This allows the missing points to 
be interpolated in the trajectory. Fig. [4] illustrates an example of 
a time series for a selected BP and its corresponding trajectory 
as computed by the tracking algorithm. Time proceeds row by 
row and from left to right. 

Normally, the standard image processing and tracking al- 
gorithms locate a BP with an uncertainty of l/N of a 
pixel (Hunter et a l] 1201 ih, where is the diameter of the BP 
in pixels. Crocke r & Grie j ((1996) did an error analysis for the 
tracking technique employed here and came to the conclusion 
that the error for locating the objects (BPs in our case) can have a 
minimum value somewhat better than 0.05 pixel. However, con- 
sidering the fact that the size and intensity variation of such BPs 
may aff'ect this precision (H unter & Weeksll201 ll) . we (over-) es- 
timate thatouralgorithi^^ to locate the BPs within half a 
pixel (ICrocker& Griejll996h . 



3.2. Properties of the BPs 



For each BP in every frame, in addition to the location, a num- 
ber of further parameters were determined, as described below. 
Some quantities, such as intensity and size (diameter), are de- 
pendent upon the threshold used to identify the BP or compute 
these parameters, since the threshold defines which pixels are in- 
cluded in the measurements. This point should be kept in mind 
when the results obtained in the following are compared with the 
other, previously found values. The diff'erence in the spatial res- 
olutions of diff'erent observations, however, is also considerable 
and very likely has a significantly stronger eff'ect. 
A BP is considered to be detected if its intensity (i.e., contrast) 
in a Ca II H image, at any given time, rises above the mean in- 
tensity of the entire frame ({Ica})- Although this criterion may 
cause the relatively low contrast BPs (i.e., the ones which are 
less bright than the mean intensity of the whole image, but have 
locally high contrast) to be missed, it reduces the chance of false 
detections. 

We discard interacting BPs to avoid falsely assigning diff'erent 
features to the same BP in the course of its lifetime. This crite- 
rion may bias the lifetime distribution to shorter lifetimes since 
the longer lived BPs are statistically more likely to interact with 
other features. In addition, limiting ourselves to non-interacting 
BPs may introduce a bias in the proper motion speeds, although 
this is not clear. Thus, such non-interacting BPs may move 
slower since the chances for interaction increases for the faster 
moving BPs that travel a greater distance. Conversely, the pro- 
cess of interaction could slow the BPs down, so that the selected 
BPs may be biased towards higher proper motions. 
The reader should keep all such biases (caused by the selection 
criteria) in mind when the properties of our selected BPs are 
compared with previous findings (where diff'erent selection cri- 
teria were applied and which therefore suff'er from diff'erent bi- 
ases). 
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3.2.1. Size 

We apply two methods to measure the BP sizes. 
Method (1): The size (diameter) of a BP is give n by the radius of 
a circl e with the same area as the BP. Following lCrocker & Griej 
(Il996l) . the effective radius R of a. Sunrise Ca ii H BP with 
arbitrary shape is calculated as the root mean square (rms) of 
intensity-weighted distances of all pixels within a BP from its 
location i.e., its centre of gravity of intensity obtained as in 
Sect. 13.1.21 This calculation is carried out within the circular 
window defined in Sect. l3.1.2l which, by construction, is always 
larger than the BP (and the BP's intensity will have dropped 
to below the threshold in all directions within this circle). 
This radius is misleadingly called "Radius of Gyration" in the 
particle tracking algorithm. It is misleading because it is not 
related to gyration or rotation of a particle or a BP. Given its 
success in following point-like sources on a noisy and variable 
background, we believe that this is an appropriate method for 
measuring the size of Ca ii H BPs since it remains largely 
unaff'ected by a sudden brightening due to, e.g., a passing wave, 
etc. Hence, it is our method of choice. 

Method (2): Sizes of BPs are estimated based on the Full Width 
at Half Maximum (FWHM) of the intensity peak at the location 
of the BP. 

The FWHM of a 2D Gaussian fit to the intensity profile of 
photospheri c BPs has been widely used for estimating thei r 
sizes (e.g.. Title & Berged ll996l: [Sdnchez Almeida et aLlfeool) . 
The 2D Gaussian functions must deal with the narrow and dark 
intergranular lanes in which the photospheric BPs typically 
reside. Unlike the photospheric BPs, the Ca ii H BPs are 
embedded in a relatively dark internetwork area with passing 
bright waves. In addition, the higher contrast of Ca II H images 
compared to photospheric passbands (e.g., G-band) as well as 
the fact that the Ca ii H data sample higher layers of the solar 
atmosphere, imply a larger size and more shape variation of the 
detected BPs. In particular, the shape variations decrease the 
quality of the fits with any specific 2D Gaussian function to the 
identified Ca II H BPs. An alternative approach is to measure 
the area across the intensity profile of the BP when it drops to 
half of its intensity peak (maximum) value. Assuming a roughly 
circular shape of the BPs and hence a circular level set for the 
Gaussian distribution, the (mean) FWHM of a BP is calculated 
as the diameter of a circle with the same area as measured. 
The latter definition (interpretation with a FWHM) results in 
diameter that are on average 21% larger and more variable 
than the ones obtained from the tracking algorithm. We prefer 
the definition based on the tracking algorithm (method 1), 
since it deals better with sudden brightenings of the image 
when travelling bright waves cross the BP. At such mo- 
ments method 2 can give an arbitrarily large size. Therefore, 
the size values used in this paper are based on the first definition. 



3.2.2. Lifetime 

The birth and death times of many BPs were not observable in 
our image sequences mainly because of the relatively short time 
series as well as the relatively small field-of-view. Therefore, the 
lifetime was computed only for those BPs that were born and 
died within the course of the time series and which did not enter 
or leave the field-of-view during the time series. A BP is con- 
sidered to be born when its brightness first rises above (Ica) and 
is considered dead when its intensity is observed to permanently 
drop below (/ca), i.e., it does not rise again to above this detect- 



ing threshold for 160 sec (i.e., 2 times our lifetime threshold) in 
the following frames. These additional frames are considered in 
order to avoid underestimating the lifetimes of BPs with repeti- 
tive brightenings, e.g., the persistent flashers. If a BP reappears 
within 30 pixels (i.e., about twice our BP's size threshold) of its 
original location within 160 sec we consider it to be the same BP. 
If these conditions are not met then we assume the new bright- 
ening to belong to a new BP. 

The relatively short time series biases the measurements towards 
shorter lifetimes and hence the lifetime distribution must be cor- 
rected (see Sect. 14.3.41) . Two possible sources of uncertainty of 
the measured lifetime are the exposure time of the SuFI Ca ii H 
images (i.e., 1 sec) and the observing cadence. The latter pa- 
rameter difl'ers for difl'erent time series used in this study and 
thus we consider an average value of the cadences presented in 
Table [U weighted by the number of BPs identified in each image 
sequence. Therefore, the Ca ii H BPs' lifetimes reported in this 
paper have an overall uncertainty of 9 sec. 

3.2.3. Intensity 

The mean value of the intensity inside the area of each detected 
BP is referred to as its absolute intensity. The intensity val- 
ues are measured from the original phase-diversity reconstructed 
Ca II H images, by referencing their accurate positions and sizes 
obtained from the restored images (described as in Sect. 13.1.11 
Sect. [3T2I and Sect. [TO) . The contrast of the BP is then deter- 
mined by normalising its intensity to the mean intensity of the 
entire Ca 11 H frame that the BP is observed in, i.e., we form the 
contrast relative to the average quiet- Sun. 



3.2.4. Horizontal Velocity 

The instantaneous horizontal velocity (v) was determined by tak- 
ing the diff'erence between the locations of a BP (known to sub- 
pixel accuracy) in two consecutive frames and dividing by the 
time between frames. The statistical uncertainty in measuring 
this frame-to-frame horizontal velocity (cr^) was computed from 
the error propagation analysis. 

We can then assess the mean value of the horizontal velocity 
(v) and its uncertainty for each B P, consisting of n time steps 
(frames), using the weighted mean (IWaU & Jenkinsll2003l) . 



<v) 



(3) 



where vj and ctj^ are the 7th determined instantaneous horizontal 
velocity (v) and its variance, respectively. The best estimate of 
the uncertainty of (v) is then. 



1 



^<v> 



(4) 



We also use Eqs. Q and (|4]) to calculate the total mean value of 
the horizontal velocity and its uncertainty. Note that we do not 
take the normal average of the velocity values to determine the 
mean velocities (for each BP or the mean over all BPs), since 
diff'erent velocity values have diff'erent uncertainties. This latter 
eff'ect arises because of diff'erent distances a BP may travel in 
equal time steps and/or diff'erent time steps due to (1) a variable 
cadences within an image sequence caused by missing or dis- 
carded images, and/or due to (2) a temporarily undetected BP. 
Such an analysis results in an average uncertainty (weighted av- 
erage over all BPs) of 0.02 km s~^ in the horizontal velocities. 
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3.3. Magnetic Field 

The magnetic properties of the Ca ii H BPs were investigated 
by comparing with maps of CP and LP determined as described 
in Sect. 12.21 In order to avoid noise introducing false detections, 
isolated single pixels lying above 3cr were not counted as sig- 
nal. The magnetic regions lying above 3cr are outlined in Fig. [4] 
by red and yellow contours for opposite polarities of the CP. A 
Ca II H BP is then considered magnetic if it overlaps with the 
computed photospheric magnetic patches for more than 50% of 
its lifetime. The magnetic field is generally still present at times 
when it is not associated with the BP but either lies below the 3cr 
threshold, or just fails to overlap with the BP. The BP shown in 
Fig. [4] seems to be almost completely inside a magnetic patch of 
negative polarity for most of the time. However, there are times 
(e.g., at the time 3041 sec from UT midnight) when this BP is 
observed just outside the boundary of the magnetic patch. 
The peak values of CP and LP over a BP's area is considered 
as corresponding polarisation degrees of the BP measured from 
the spatially smoothed non-reconstructed data (obtained as de- 
scribed in Sect. 12.21) . 

3.4. Inversions 

We use the results of SIR code (Stokes Inversion ba sed on 
Response functions : lRuiz Cobo & Del Toro IniesTalll992h inver- 
sions of the IMaX Fe I 5250.2 A Stokes vector carried out by 
L. R. Bellot Rubio. The code numerically solves the radiative 
transfer equation (RTF) along the line of sight (LOS) for the 
Zeeman-induced polarisation of light under the assumption of 
local thermodynamic equilibrium (LTF). The code deals with 
all four Stokes parameters (/, Q, U, V), any combinations of 
which are fitted for any arbitrary number of spectral lines. 
Then, the diff'erences between the observed and synthetic Stokes 
profiles are iteratively minimized by modifying an initial model 
atmosphere. 

Details of the employed SIR inv ersions can be found in the 
paper by lGuglielmino et aP (l201lh . 

The Harvard Smithsonia n Reference Atmosphere 
(HSRA; iGingerich et all Il971>) is used as the initial model 
atmosphere. The temperature has two nodes and the inver- 
sions recover the temperature stratification in layers between 
logir) - 4 and (where r is the optical depth of the continuum 
at 5000 A). Stray light is not taken into account for this inversion 
and a magnetic filling factor of equal to unity is assumed. This 
implies that the magnetic field strength from the SIR code is a 
lower limit of the true value. 

4. Results and Discussion 

A total of 107 Ca ii H BPs were identified in 6 diff'erent im- 
age sequences obtained by the SuFI instrument, where each BP 
is counted only once during its lifetime. We found a relatively 
low BP average number density of 0.03 (Mm)~^ corresponding 
to ^ 1.8 X 10^ BPs at any given time on the whole solar sur- 
face. This number is roughly an order of magnitude smaller than 
the number of granules on the solar surface at any given time. 
This low number can be due to the rather restrictive identifica- 
tion criteria, in particular that we consider only small isolated 
BPs (comparable in size with the Sunrise spatial resolution) that 
do not merge or split while they migrate. These criteria restrict 
our sample mainly to internetwork BPs. The parameters inten- 
sity contrast, horizontal velocity and size, were calculated for 
each BP at every time step. The lifetime was measured only for 



47 BPs whose birth and death times were both observable. The 
BP is considered magnetic if the average of the local Stokes V 
signal within a 3 x 3 pixel box around its corresponding location 
is higher than the 3cr level over at least 50% of the detected BP 
lifetime. 



4.1. A Case Study: Tracking an individual Ca ii H BP 

Here, we study the evolution of one of the identified BPs in de- 
tail. Snapshots of the Ca ii H intensity in the region around the 
BP as well as overlaid contours of CP are plotted in Fig.|4l The 
BP moves within an area of roughly 1x2 arcsec^ in the course of 
its lifetime of 1321 ± 13 sec (see its trajectory in the bottom right 
frame). The BP is located almost always inside a magnetic patch, 
except for a few moments such as the one that occurs at time 
3041 sec (all times are given relative to UT midnight). At this 
time, the BP lies mainly outside the boundary of the magnetic 
patch, however, a small part of its area overlaps with CP > 3 
signal. Note the irregular intervals between the frames. Fig. [5] 
illustrates the temporal behaviour of physical parameters of the 
BP. It has a mean intensity of 1.65 ± 0.01 {Ica) and moves with 
a mean proper motion velocity of 1.86 ± 0.08 km s~^. 
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Fig. 4. Temporal evolution of a selected BP (Lifetime: 1321 sec), 
observed on 9 June 2009 in the Sunrise/SuFI Ca ii H passband. 
Time proceeds row by row and from left to right within each 
row. The red and yellow contours outline magnetic regions with 
opposite polarities of the line-of-sight magnetic field component 
and match the 3cr level in Stokes V. In the bottom-right panel, 
the trajectory of the BP is plotted. The arrow marks the sense of 
the motion. The time from UT midnight, in seconds, is shown in 
the upper-right corner of each frame. Note the irregular intervals 
between frames. 
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The frame-by-frame intensity variation (thin curve) is plotted in 
panel (a). The intensity curve smoothed by a boxcar averaging 
of 5 frames is overplotted (red in the online version). In order to 
estimate the effect of transient brightenings, e.g., passing waves, 
on the bp's intensity variation, we overplot the variation of the 
averaged intensity within the BP (smoothed by a 5 -frames box- 
car averaging; blue dashed curve) measured in the restored im- 
ages, i.e., the images from which the larger wave-like features 
are subtracted (see Sect l3.1.l1) : the curve is in arbitrary units. A 
comparison between the two smoothed intensity curves validates 
the intrinsic intensity variation of the BP. Panel (b) represents 
the variation of horizontal velocity smoothed over 5 frames by 
weighted averaging, described in Sect. 13.2.41 (thick curve; red in 
the online version). The thin vertical solid lines indicate the er- 
ror bars. The values of intensity and horizontal velocity vary by 
a factor of nearly 2 and 3, respectively. 

Close inspection of panels (a) and (b) reveals an anti-correlation 
between these two quantities for the BP. The vertical dashed 
lines mark the times when the intensity maxima and minima fall 
together with minima and maxima in horizontal velocity, respec- 
tively. However, two peaks in the horizontal motion plot and the 
corresponding dips in intensity plot, indicated by vertical dot- 
dashed lines, show a small time delay. An anti-correlation be- 
tween the intensity and horizontal velocity is also observed in 
some other BPs, sometimes also with a small time delay. All 
BPs do not show it. 

Panel (c) in Fig. [5] illustrates the evolution of the BP's diameter, 
which is found to be constant to within 10 - 15 %. 
The CP (obtained as described in Sect. 12.21 and normalized to 
the continuum intensity (7^) averaged over the full IMaX FOV), 
plotted in panel (d), is always above the 3cr noise level, which 
demonstrates the magnetic origin of this BP. 
The variation of total linear polarisation (LP; see Sect. 12. 2K nor- 
malised to (/c), is plotted in panel (e). The dashed and dotted 
lines in panel (e) show the 3cr and 2cr noise levels. The LP is 
almost always below the Icr level (indicated by the dot-dashed 
line), implying that the linear polarisation signal detected at 
the position of the small BP is not significant at the Icr level. 
Nevertheless, inversions return a strongly inclined field for such 
a BP with inclination angles between 78 and 90° in the course of 
its lifetime. This result, however, is probably aff'ected by noise 
as we find by comparing with inclinations derived from the lo- 
cations of the BP in images sampling diff'erent heights in the 
atmosphere (Jafarzadeh et al., in preparation). 
The values of the field strength multiplied by the absolute value 
of the cosine of the inclination angle \B cos(7)| (i.e., the unsigned 
longitudinal component of the magnetic field) returned by the 
SIR inversion code are plotted in panel (/). We plot the varia- 
tions of \B c os(7)| since this value is retrieved more reliably than 
B (see Sect. 14.3.51) . \Bcos(y)\ varies between 108 and 361 G as- 
suming that the magnetic feature is spatially re solved. If this par- 
ticular BP is spatially unresolved by Sunrise (iRiethmiiller et al.l 
2012), which is likely given the low CP value, these field 
strength values are the lower limits. Fig. [6l shows a \Bcos(y)\ 
map of the BP at time 3304 sec from UT midnight. \Bcos(y)\ 
at the position of this small internetwork BP reaches a maxi- 
mum value of 322 G at this point in time. The over-plotted solid 
contours show the field strength levels of 100, 200 and 300 G. 
The location of the Ca ii H BP is overlaid (dashed contour). 
The maximum intensity of the the BP, indicated by a cross, is 
off'set by 0.05 arcsec (20% of the diameter of the BP) relative 
to the location of the highest \Bco^{y)\. This off'set may repre- 
sent a small inclination of the flux tube from the vertical direc- 
tion, since the Ca ii H brightness and \Bco^{y)\ are determined 
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Fig. 5. Evolution of {a) the intensity contrast, (b) horizontal ve- 
locity, (c) diameter, {d) CP, a measure of the strength of Stokes 
y, {e) total net linear polarisation LP, and (/) field strength of the 
selected BP shown in Figs. [T] and IH CP and LP are normalised 
to Stokes / continuum. In panel (a) and (/), the thin curves re- 
spectively represent the time variation of average intensity and 
the unsigned longitudinal magnetic field within the BP while the 
smoothed thick lines (red in online version) show a 5 -frames 
boxcar average. They are both normalised to the mean intensity 
of the image. The (blue) dashed curve in panel {a) represents 
the smoothed average intensity within the BP, in arbitrary unit, 
from the restored images (see main text). Evolution of horizon- 
tal velocity, smoothed within 5 frames by weighted averaging, is 
plotted as the thick line (red in the online version) in panel (b) 
while the thin vertical solid lines indicate the error bars. The ver- 
tical dashed lines in the two upper panels serve to guide the eye 
to the general anti-correlation between intensity and horizontal 
velocity (see main text). The 3cr noise levels for both CP and LP, 
are drawn as dashed-lines in panels (J) and {e), while the dotted 
and dot-dashed lines in panel {e) are the 2cr and Icr noise level^ 
respectively. The unsigned longitudinal magnetic field depicted 
in panel (/) was computed from the SIR inversions (see main 
text). 
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at different heights. We will consider inclinations of magnetic 
elements in a separate paper. A small observing time difference 
between the SuFI Ca ii H and IMaX images can also be a reason 
for such an offset. 



■ I I I I I I I I I I I I I I I I I I I I I I I I I I 



1.0 - 
0.8 - 



o 

^ 0.6 h 



0.4 - 
0.2 - 





I I I ■ ■ I I I ■ ■ 1— ■ ■ ■ I I I I I I I I I Jb 

0.2 0.4 0.6 0.8 1.0 1.2 



|— 1322 



-249 



- 176 



130 



30 



O 

o 
o 

CO 



[arcsec] 



Fig. 6. Map of the unsigned longitudinal component of mag- 
netic field, \Bcos(y)\, centred on the BP shown in Fig.|4]at time 
3304 sec from UT midnight. The solid contours indicate the 
height-independent \B cos(r)| at 100, 200 and 300 G. The Ca ii H 
BP 's area is indicated by the dashed contour. 
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Fig. 7. T emporal evolution o f a BP similar to the persistent 
flasher of iBrandt et al.l (Il992h . Time proceeds row by row from 
left to right in each row. The yellow and red contours indicate 
magnetic patches of opposite polarities. The time in seconds, 
from UT midnight, is given in the upper-left corner of each 
frame. 

The flashers differ from the normal BPs mainly in that they (the 
flashers) tend to be less bright. The persistent flashers are on 
average 47% less bright than the other BPs studied here (see 
Sect. 14.3.1b . This causes the occurrence of intensity variations 
(^^ of a factor of 3) around our detection limit and hence the 
intermittent appearance of the BP. We conclude that persistent 
flashers are less bright, but otherwise normal magnetic BPs. 



4.2. Persistent Flasher 

Among the 107 a nalysed BPs 7 are si milar to those named per- 
sistent flashers by lBrandt et al.l (Il992h . Their main difference to 
the majority of the selected BPs is the occurrence of intensity 
variations about the detection limit ({Ica}) which give rise to the 
apparent succession of flashes. 

An example of such a BP is shown in Fig. [71 which shows the 
Ca II H intensity at 6 times within an area of 2 x 2 arcsec^. The 
contours outline magnetic regions with CP above the 3cr level. 
The temporal variations of different parameters of the BP are 
plotted in Fig. [HI From panel (a) it is evident that between two 
bright phases the intensity of the BP drops below the intensity 
limit of the detection giving it the appearance of flashing, al- 
though its intensity variations are no larger, in a relative sense, 
than those of the BP displayed in Fig [3 Interestingly, an anti- 
correlation between the intensity and the horizontal velocity in 
panel (b) is also seen for this persistent flasher. The flasher trav- 
els with a mean proper motion velocity of 2.82 ± 0.06 km s~^. 
The size is constant within a range of 10% (panel (c)). We can- 
not rule out that the large excursion of the size at the beginning 
of the plotted time series is due to a false detection at the very 
beginning of the BP's identification. Similar to the other BPs un- 
der study, CP always lies above the 3cr noise level, also for the 
persistent flasher; see Fig. [8]i. Again, the LP signal lies almost 
always below the cr level (dot-dashed line). Therefore, the same 
conclusion as for the BP studied in Sect. 14.11 can be made here 
of the field inclination. The \Bcos(y)\, plotted in panel (/), has 
a mean value of 190 G, smaller than the majority of the normal, 
non-flasher BPs under study. 

In this investigation, the few "persistent flashers" (6% of all iden- 
tified BPs) clearly have a magnetic origin, with horizontal veloc- 
ities of 2.4 km on average. This makes the persistent flashers 
about equally mobile as the other Ca ii H BPs (see Sect. 14.3.21) . 



4.3. Statistical Studies 

Table [2l summarises the number and mean values of the BPs 
studied here. The mean lifetime value is to our knowledge the 
longest obtained for Ca ii H BPs so far (see Sect. 14. 3. 41) , although 
it is a lower limit due to the finite lengths of the employed time 
series. This reflects the excellent quality of the Sunrise/SuFI ob- 
servations. The mean diameter is influenced by our criteria to 
only concentrate on the isolated and point-like features which 
do not display any internal fine structure. 

Table 2. Summary of averaged properties of the Ca ii H BPs 



Parameter 


Mean value ^ 


No. of BPs 


Intensity 


1.48 + 0.3 (Ica) 


107 


Horizontal Velocity 


2.2 ± l.^kms-^ 


107 


Diameter 


0.2 + 0.02 arcsec 


107 


Lifetime 


673 + 526 sec 


47 


CP^ 


0.32 + 0.21 % 


53 


|5cos(7)|^ 


142 + 87 G 


53 



^ The uncertainty estimates were calculated from the standard devi- 
ations of the parameters' distributions. For formal errors see main 
text. 

^ From SuNRisE/IMaX 

CP: Circular polarisation degree, described as in Sect. I2.2[ 
|5cos(7)|: unsigned longitudinal component of magnetic field. 



Distributions of the intensity, horizontal velocity and size of 
all 107 BPs, as well as the lifetime of the 47 BPs for which it 
is available, are presented in Fig. [9l Before discussing these 
properties in detail, we point out that all of the BPs for which 
simultaneous polarimetric observations were available (^ 50% 
of all BPs) are clearly associated with magnetic fields (see 
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of the common origin and the similarity in appearance of 
the magnetic BPs in G-band, TiO-band and Ca ii H/K pass- 
bands ( de Wijn et al. 2008), we have also included G-band and 
TiO passbands, which normally sample lower regions of the 
solar atmosphere than Ca ii H/K. However, the sampled heights 
depend on the filter widths. Thus, the measured Ca ii H/K signal 
also gets contributions from the lower layers of the atmosphere. 
The average horizontal speed of non-magnetic K-grains is also 
listed in this table in order to illustrate the diff'erence to magnetic 
BPs. 




1.0 1.5 2.0 
Intensity/I^e 



2 4 6 8 10 12 
Horizontal Velocity [km/s] 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 




(C) 





\ 


















\ 




























S- - 









100 140 180 
Diameter [km] 



200 600 1000 1400 
Lifetime [s] 



Fig. 9. Statistical properties of the detected Sunrise Ca ii H BPs. 
(a)\ Distribution of mean intensity contrast for all 5677 individ- 
ual measurements (black, outlined histogram) and of the time 
averaged contrast of all 107 BPs (grey, filled histogram). (b)\ 
Distribution of all 5570 individual measurements of horizontal 
velocity (black, outlined histogram) and of the values averaged 
over their respective lifetimes of the 107 studied BPs (grey, filled 
histogram), (c): Histogram of mean size (diameter) of all 107 
BPs (grey, filled histogram) and all 5677 individual measure- 
ments (black, outlined histogram). The vertical dashed line in- 
dicates the spatial resolution achieved by the Sunrise telescope. 
(d): Histogram of lifetimes of 47 BPs (black, outlined histogram) 
and of corrected lifetimes (grey, filled histogram; see main text). 
The vertical solid lines indicates the mean values of the distribu- 
tions. The fits (dot-dashed and dashed curves) are explained in 
the main text. 



Sect. 14.3.5b . It is useful to bear this in mind when reading the 
following sections. 

The diff'erence in the spatial resolutions and scattered light 
between diff'erent sets of data in the literature is considerable 
and may well be a dominant cause of discrepancies between 
various studies. In addition, observing in diff'erent wavelengths 
and hence sampling diff'erent heights in the solar atmosphere 
affects the results. The same is true for the variation in the 
lengths of observed image sequences. Finally, diff'erences in 
definition of BPs may also influence the results (e.g., thresholds 
used). 

IVLean values of BP physical parameters determined here are 
compared with past findings in Table O In this table, because 



4.3.1. Intensity 

The distribution of all 5677 measured instantaneous intensity 
contrasts of the 107 BPs is represented by the outlined his- 
togram in Fig. [9]2. Its exponential fit with an e-folding width of 
0.22 (Ica) is overlaid as a red, dot-dashed line (fit limited to the 
range of 1.45 to 2.7 (Ica))- 

The intensity contrast distribution shows a lower limit of 
0.69 (Ica) and the brightest detected BP reaches a maximum 
value of 2.86 (Ica)- The mean contrast value, averaged over all 
BPs and all time steps is 1.48 (Ica)- The few intensity values 
less than (Ica) belong to persistent flashers whose intensity 
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drops below the detection limit ((Ica)) at some moments of 
time, so that the average is low, although in between they are 
clearly visible. 

The filled histogram in Fig.[9]2 is the distribution of the lifetime 
average of intensity contrast of each individual BP of all 107 
BPs. 

A comparison with values in the literature is summarised 
in Table [3] Such a comparison is not straightforward because 
observations at diff'erent wavelengths, with diff'erent filter widths 
and diff'erent spatial resolution give a range of contrasts. Also, 
none of the non-SuNRiSE studies of Ca ii H or K BPs provides 
their average contra sts. 

iBerger et all (Il995h showed that the intensity of their GBPs, 
on average, reaches 1.27 relative to the quiet-Sun average. A 
ra nge of 0.8 - 1.8 (Ips) ^yas repo rted for the same wavelength 
by ISanchez Almeida et al.l (l2004l) at an estimated spatial reso- 
lution of 135 km on the Sun. Note that the lower boundary of 
this range means that these BPs are darker than (Iqs ) implying 
that the GBPs were bright only with respect to the intergranular 
lanes surrounding them. A mean value 1.17 (Iqs) was obtained 
bviMostl et al. ( 2006) for magnetic GBPs. 
iRiethmiiller et al.l ([2010) studied the peak contrast of the BPs 
observed in diff'erent passbands of the Sunrise observatory 
between 2140 A and 5250.2 A, including Ca ii H 3968 A. They 
found mean values of the peak contrast equal to 1.11, 1.35, 1.52, 
1.60, 1.89 and 2.31 (Iqs) for 5250, 3120, 3000, 3880, 3968 and 
2140 A, respectively. Their mean contrast value of the Ca ii H 
BPs is somewhat larger than the 1.48 we measured in this study. 
This diff'erence partly reflects the difl'erent criteria employed to 
identify the BP studied (we restrict ourselves to smaller BPs). 
Moreover, they measured the peak contrast of each BP, which 
is expected to be higher than the contrast averaged over a BP's 
area as we determined. 



4.3.2. Horizontal Velocity 

The distribution of the instantaneous proper motion velocity of 
all individual 5570 frame-to-frame measurements made in this 
study is given by the outlined histogram in Fig. (2?. The red, dot- 
dashed line is the exponential fit to this histogram (for velocities 
larger than 1.8 km s~^) with an e-folding width corresponding 
to a horizontal velocity of 1.8 ^"^ However, lower velocities 
(< I km s~^) are under-represented relative to the exponential fit, 
suggesting that the BPs tend not to stay at a given spot. The dis- 
tribution displays a tail reaching up to 10 ^"^ There are a few 
measurements at higher velocities which lie outside the range of 
this histogram reaching upto 15.5 km s~^. Using the values of 
temper ature, density and gas pressure of the VAL-C atmospheric 
model ( IVernazza et al.ll 198 11) at the average formation height of 
the Sunrise Ca ii H BPs (500 km), the sound speed is estimated 
to be between 6.5 and 7 km . Thus, Fig. [2? implies that about 
3.5% of the time the BPs' motion is supersonic. 
The distribution of the mean value (lifetime average) of the hori- 
zontal velocity of each BP is over-plotted (filled histogram). The 
exponential fit to this histogram (with an e-folding width corre- 
sponding to a horizontal velocity of 1.3 ^"^) is plotted as a 
blue, dashed line. Power law fits were found to be unacceptable 
for both these distributions. The statistics reveal that the mean 
horizontal velocity of the BPs studied in this work varies be- 
tween 1.2 and 6.6 km (mean values for individual BPs) with 
an average value of 2.2 km (averaged over all BPs). 
We remind the reader that the value of the frame-to-frame ve- 



locity is biased by the measurement error in positioning the BPs 
(see Sect. 13. 1.2b . An upper limit on the velocity induced by this 
error is about 0.7 km s~^. The eff'ect of this error is greatly de- 
creased when computing the velocities averaged over the life- 
times of individual BPs, as mentioned above (see Sect. l3.2.4b . 
As summarised in Table [3l the measured horizontal veloc- 
ities of the Ca ii H BPs are larger than the velocities 
of magnetic BPs in the photosphere, wh ich are quoted as 



0.06 km s'^ (Ivan Ballegooiienet all Il998l) . 1 - 1.4 km s 



(iMostl et al.ll2006h . although higher horizontal velocities of up 
to 5 km for a few indicate d BPs have b een reported in the lit- 
erature (iBerger & Titld[T99 6). Note that van Ballegooijen et"aD 
(19981) determined horizontal velocities by tracking corks so 
that these refer to averaged motions of clusters of BPs instead 
of isolated ones. However, our va lue is significantly sm aller 
than the ave rage of 6.6 km s -^ (IWellstein etal.1 1 1998h and 
49 km ^"^ ( Steifens et al.lll996l) reported for the Ca ii H inter- 
network K-grains as well as 7 - 10 ^"^ horizonta l velocities 
of the network BPs found bv IWeUstein et al.l (Il998l) . Whereas 
the somewhat larger values we find compared to the GBPs may 
be acceptable (e.g., caused by swaying of magnetic flux tubes), 
the huge range of values quoted for the chromosphere (Ca ii H 
and K observations) shows how difficult it is to compare these 
values. Obviously, the diff'erent measurements refer to very dif- 
ferent kinds of features. In particular, K-grains refer to a wave 
phenomenon and cannot be compared with the magnetic fea- 
tures studied here. There still remains the discrep ancy with th e 
large value for network BPs published by Wellstein et al.l (1 19981) . 
Possible explanations for this are: (a) Network BPs (NBPs) 
move much faster than internetwork BPs (IBPs). (b) Because 
Wellstein et al. (1998) had no magnetic information, their mea- 
surements may have been aff'ected by misidentification of pass- 
ing waves even in the case of NBPs. We obtained stable results 
only after applying a careful filtering technique, (c) Proper mo- 
tion velocities derived from ground-based observations tend to 
be exaggerated due to diff'erential seeing-induced deformations 
in diff'erent frames of the time series. 

Since we are unaware of any studies (e.g., in the G-band) sug- 
gesting higher horizontal velocities of network BPs than inter- 
network BPs, and due to the m uch larger horizontal velocities 
found by IWellstein et al.l (1 19981) than of photospheric BPs, we 
conclude that explanations (b) and (c) are more likely than (a). 



4.3.3. Size 

The distribution of the BP diameters obtained from all 5677 indi- 
vidual measurements are represented by the outlined histogram 
in Fig. [9]:. The BPs' diameters range between 85 and 195 km or 
^ 0.11 and ^ 0.27 arcsec, with a mean value of 150 ± 15 km 
(^ 0.2 ± 0.02 arcsec). 

The filled gray histogram represents lifetime averages. Both his- 
tograms exhibit a rather narrow peak centred at 150 cor- 
responding to a diameter of ^ 0.2 arcsec. This narrowness of 
the peak reflects our size criterion in which a Sunrise Ca ii H 
BP was defined as a bright, point-like feature with a maximum 
diameter of 0.3 arcsec at any given time. On the other hand, the 
BPs' size distribution extends down to the resolution reached by 
Sunrise of roughly 100 km in the best time series. Therefore, 
the typical diameter of the BPs studied here is between 0.15 and 
0.25 arcsec. Possible explanations for the missing smaller BPs 
in our statistics (down to the theoretical resolution limit) are: 
(a) Resolution of the Ca ii H images is somewhat larger (e.g., 
0.2 arcsec) than the best resolution reached by Sunrise, (b) The 
smaller features are not bright enough to be identified by our 
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Table 3. Comparison of the mean values of properties of BPs obtained in this study with values in the literature. 



Reference 


Telescope 
(Diffraction 
Limit) 


Passband 
(FWHM) 


Pol. 
Info." 


Feature*^ 


Intensity 

KIqs)] 


H. Velocity 

[km s~^] 


Diameter 

[arcsec] 


Lifetime 
[^] 


This study 


Sunrise (Of! 10) 


CaiiH(1.8 A) 


Yes 


IBP 


1.48 


2.2 


0.2 


673 


^iethmuller et al.^ (2012) 


Sunrise (a.'lO) 


CN (8.0 A)'^ 


Yes 


MBP 


_ 


_ 


0.45 


_ 


Kevs et al. (2011) 


DST (Of! 12) 


G-band (9.2 A) 


No 


MBP 


_ 


1.0 


_ 


91 


Abramenko et al.^ (2010) 


NST (QfUl) 


TiO (10 A) 


No 


IBP 


_ 


_ 


^0.16 


< 120 


Q-ockett et al (2010) 


DST {Of: 12) 


G-band (9.2 A) 


No 


MBP 


_ 


_ 


0.31 


_ 


Riethmuller et al. (2010) 


Sunrise (a^O) 


Ca iiH(1.8 A)'^ 


Yes 


MBP 


1.89 


_ 


_ 


_ 


Mosti et al. (2006) 


SST (Of! I) 


G-band (11.6 A) 


No 


MBP 


1.17 


1.11 


0.28 


263 


de Wiin et al. (2005) 


DOT (a'2) 


Ca II H (1.3 A) 


No 


IBP 


_ 


_ 


_ 


258 


Sanchez Almeida et al. (2004) 


SST (a'l) 


G-band (10.8 A) 


No 


GBP 


0.8-1.8 


_ 


0.25 


225 


Nisenson et al. (2003) 


DOT (a'2) 


G-band (10 A) 


No 


GBP 


- 


0.89^ 


- 


552 


Berger& Title (2001) 


SVST (a'2) 


G-band (12 A) 


Yes 


GBP 




1-5^ 


0.145 




Wellstein et al. (1998) 


VTT (a'3) 


Ca II K (0.3 A) 


No 


NBP 




6.6 






Steffens et al. (1996) 


VTT (a'3) 


Ca II K (0.3 A) 


No 


K-grain 




50 






Berseretal. (1995) 


SVST (a'2) 


G-band (12 A) 


No 


GBP 


1.27 




0.35 




Muller et al., (1994) 


Pic du Midi (a.'25) 


5750 A (60 A) 


No 


NBP 




1.4 




1080^ 


Soltau (1993) 


VTT (a'3) 


Ca II K (0.3 A) 


No 


IBP 




2.5 






Muller (1983) 


Pic du Midi (a.'25) 


5750 A (60 A)'^ 


No 


IBP 








540 



" Polarisation (magnetic field) information 

^ IBP: Internetwork Bright Point; NBP: Network BP; MBP: Magnetic BP; GBP: G-band BP 

^ See main text for other passbands 

^ rms velocity 

^ From Berger & Titld (119961) 

^ From|Mullerj (Il983h , who found this value for NBPs. 



detection procedure, (c) The smaller features are not long-lived 
enough to be considered as a BP based on our restricted criteria 
(see Sect. [3]). 

Comparing the diameters obtained here with those reported in 
the literature is not straight forward, since they depend on the 
spatial resolution of the employed observations, the level of so- 
lar activity in the observed area (internetwork, network, plage, 
etc.) and on the technique used to determine the size. In this 
study, in order to avoid the effect of fine structure, we have set 
an upper limit on the size, so that we do not include the larger 
features that are also present in the Sunrise data. We review the 
reported values of the BP sizes observed in the photosphere and 
the low chromosphere in Table [3l However, it is important to re- 
member that the magnetic field of magnetic elements expands 
with height. Therefore, under the assumption that the brightness 
structures scale with the size of the magnetic features, a BP seen 
in the Ca ii H passband should be larger than in the G-band, and 
larger still than in the TiO band (which corresponds essentially 
to the continuum outside sunspot umbrae) under the same ob- 
serving conditions. 

Using the network flux-t ube model of | Solankil (Il986l) and the 
thin-tube approximation (iDefouwl \l9m we estimate the ra- 
tios of the flux tube diameter at the diflTerent heights of the at- 
mosphere sampled by the red continuum (TiO band), the G- 
band (roughly 50 km higher; Carlsson et al. 20 04) and Ca ii H 
440 km higher; see Sect.O)- We obtain 1 : 1.1 : 2.6, re- 
spectively, assuming the approximate mean formation heights. 
However, these ratios represent a lower limit of the size ratios, 
since the Ca ii H BPs form in higher layers due to the eflTect 
of higher }l2v/^2R emission peaks around the Ca ii H line-core 
(see Sect. 12.11) . One should keep these expansion rates in mind 
when comparing the size of BPs in diflTerent layers. Note that the 
Wilson depression afiTects all layers and does not, to first order. 



aflTect this ratio 

Berger et^ (Il995h reported a range of FWHM between 0.17 
and 0.69 arcsec, with an average value of 0.35 arcsec for pho- 
tospheric G-band BPs. U sing the same telescope and the same 
passband, iBerger & Titld fZOOl) noted a much smaller mean 
value of 0.145 arcsec. This latter value is well bellow the the- 
oretical diflfraction limit, 0.2 arcsec, of the telescope they used 
(SVST; Swedish Vacuum Solar Telescope) and hence smaller 
than the spatial resolution of their observations. This diameter 
probably reflects the threshold they used to obtain the size. In 
any case it suggests that they probably did not resolve most o f the 
BPs in their sample. Later, ISanchez Almeida et al.l (|2004|) ob- 
tained a diameter of 0.25 arcsec on average (for GBPs), in data 
from a telescope with a higher diflTraction-limit. ^ Crockett et al.l 
(2010) employed one-dimensional intensity profiling to detect 
magnetic GBPs and found a diameter distributi on peaked at 
0.31 arcsec which is in agreement with that of iBerger et al.l 
( 199 5|). Their largest G BPs were 1 arcsec in diameter. Recently, 
Riethmuller et al. ( 2012) studied the photospheric BPs from the 
Sunrise/SuFI CN 3880 A passband and found an average diam- 
eter of 0.45 arcsec with a distribution of sizes ranging from 0.24 
to 0.78 arcsec. A relatively small mean value of 0.16 arcsec was 
reported for internetwork BPs observed in high spatial resolution 
images taken in the TiO band ( Abramenko et al.ll2010l) . An ex- 
cessively large value of the diameter has been reported for the 
magnetic Ca ii K BPs, observed in internetwork areas (ISoltaul 
1993). He reported a diameter of 2.5 arcsec and considered it 
as an indication of a preferred size for the magnetic features. 
Obviously, either he considered very diflTerent features than we 
did or his spatial resolution was very low. 
Based on the estimated expansion rates of a thin flux tube in 
the lower solar atmosphere, a photospheric GBP, i.e., the cross- 
sectional area of the flux tube in the height sampled by G-band, 
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with a diameter of 0.35 arcsec (i.e., the largest mean value 
of a GBP's diameter reviewed here) expands to a diameter of 
0.8 - 1.3 arcsec at the height sampled by Ca ii H passband, as- 
suming nearly vertical flux tubes. Such an expansion with height 
may lead to C a ii H/K BPs with diameters as large as 2.5 arcsec 
(ISoltaull 19931) corresponding t o relatively large GBPs, which, al- 
thoug h rare, are observed (e.g. jiBerger et al.l 19951: ICrockett et al.l 
l2Qll . 

4.3.4. Lifetime 

Histogram of the lifetime for the 47 BPs for which both birth and 
death times were observable are shown in Fig.[9ji. The red, dot- 
dashed line is an exponential fit with e-folding time of 526 sec, 
roughly 9 min. The condition that only BPs with lifetimes longer 
than 80 sec are considered in this study, in principle limits the 
left side of the distribution. As mentioned in Sect.O this lifetime 
threshold of 80 sec was introduced to exclude point-like bright- 
enings caused by superpos itions of waves and reversed granu- 
lation (Ide Wijn et al.ll2005l) . However, the shortest lifetime we 
found is 167 sec, so that in practice this limit should not influ- 
ence the deduced lifetimes. This result suggests that either there 
is a lower limit to the lifetime of magnetic BPs, or one of the 
other criteria by which we isolate our BPs also discarded shorter 
lived BPs. 

The longest lived BP in our sample lasted 1507 sec. This is only 
slightly shorter than the longest image sequence we have used in 
this study and therefore, the length of our image series restricts 
the longest lifetime we can determine. This causes the number 
of longer lived BP s to b e underestimated. Therefore, follow- 
ing iDanilovic et al.l (12010) . we correct the lifetime distribution 
by multiplying it with a weight of 



where m is the number of frames that the BP lives and n is the 
total number of frames in the observed time series (shown in 
shaded grey in Fig.|9]i). This results in a somewhat higher mean 
lifetime of 673 sec, i.e., roughly 1 1 min . 

In an earlier study, de Wijn et al] (l2005l) discussed that tracking 
the Ca II H BPs in internetwork regions is harder than in the net- 
work due to the larger, more dynamic granules in the internet- 
work regions that cause crashing into the flux tubes. They con- 
cluded that this interaction can disturb the processes that cause 
the brightness of BPs and hence make them invisible. They re- 
ported a mean lifetime of 258 sec for their internetwork Ca ii H 
BPs. The mean lifetime of 673 sec for internetwork Ca ii H BPs 
that we have obtained from the Sunrise data with a higher reso- 
lution, is a factor of 2.6 larger than their value and is indicated 
by a vertical line in the lifetime histogram illustrated in Fig.|9]i. 
At least a part of this difl^eren ce could be d ue to the efl'ect of 
variable seeing on the results of de Wij n et al.l (120051) . which can 
lead to significant underestimates of the lifetime. 
lAbramenko et al.l (l2010b showed that the majority (^ 98.6%) 
of photospheric BPs, identified in high resolution images, have 
a lifetime less than 120 sec. This implies either a fundamen- 
tal diff'erence between photospheric and chromo spheric BPs. 
Alternatively, it is due to the difficulty of reliably identifying in- 
ternetwork BPs in the visible continuum due to the generally low 
contrasts, particularly in the absence of polarisation information. 
Alternati vely, it may be the eff^ect of variable seeing in the obser- 
vations of lAbramenko et al.l(l2010l) . lNisenson et al found 
an averaged lifetime of 9.2 min for their isolated, network GBP, 



which is com parable with our find ing for isolated internetwork 
Ca II H BPs. de Wijn et al." ('2008') reported a mean Hfetime of 
10 min for the magnetic internetwork elements they studied in 
the quiet- sun. Although this is in a good agreement with what 
we obtained here, only 20% of their magnetic internetwork el- 
ements are related to BPs. They do not say if these BPs diff'er 
from the remaining magnetic features. In a lower resolution ob- 
servation 0.3 arcsec), Muller (1983) found a mean value of 
1080 sec for network BPs and an average lifetime of 540 sec 
for internetwork BPs, both observed in the photospheric 5750 A 
passband. In addition, they also found a range of 1 80 - 600 sec 
for the lifetime of Ca ii K BPs. Although at lower spatial reso- 
lution and aff'ected by variable seeing, the result for the internet- 
work BPs is in agreement with our findings. Note that they also 
restricted themselves to the study of "point-like" BPs by setting 
an upper limit of 0.5 arcsec for the BPs' diameters. 
Except for the res ults reported by IMulleri (119831) and 
iNisenson et al.l (l2003h . the lifetimes we find for internetwork 
BPs are significantly larger than those in the literature, both in 
the photosphere and the chromosphere. This may well be due to 
the stable conditions provided by Sunrise. An alternative expla- 
nation for the longer lifetimes obtained in the Sunrise data could 
be that through our criteria that allow a BP to disappear and ap- 
pear again, we are sometimes falsely assigning two or more fea- 
tures to the same BP in the course of its lifetime. However, it 
has only a small eff'ect on our lifetime distribution, as we found 
by recalculating the distribution without applying this criterion 
(only the few persistent flashers are afl'ected). To test this situa- 
tion, we consider the magnetic polarity of each BP for which this 
information is available. We find that every BP always maintains 
the same polarity during its lifetime. This supports our identifi- 
cations. Our criterion to restrict our BPs in size may conceiv- 
ably lead to a bias in lifetimes (if smaller BPs are longe r lived). 
The agreement between our result and the one found bv lMuUed 
(1983), who also set an upper limit on BP sizes, may support 
such an interpretation. 



4.3.5. Circular Polarisation 



Plotted in Fig. [TOl is the distribution of all (unsigned) CP val- 
ues (described as in Sect. 12.21) measured in individual snapshots 
of all 53 BPs for which simultaneous IMaX observations were 
available (see Table [B. The 3(T noise level is indicated by the 
vertical dashed line. The exponential fit to the distribution (re- 
stricted to CP > 0.2%) is overlaid (red dot-dashed line). It has 
an e-folding width of 0.21 CP [%]. The vertical solid-line indi- 
cates the mean value of CP, which corresponds to 0.32%. 
Our mean value of CP, from the spatially smoothed non- 
reconstructed IMaXdata, is smaller, by a factor of about 8, than 
that of Riethr nuTler et al.l (|2012), who measured this quantity at 
the position of their BPs from the phase-diversity reconstructed 
IMaX data. Hence, their larger photospheric CN BPs with di- 
ameters between 0.24 - 0.78 arcsec show much stronger polar- 
isation signals compared to our small-scale, low chromospheric 
elements with diameters of 0.1 1 - 0.27 arcsec. Note that due to 
the flux tube's expansion with height, our small Ca ii H BPs are 
considered to have an even smaller size at the height sampled by 
the photospheric CN passband. Our measured CP values (at the 
position of small BPs) from non-reconstructed data are roughly a 
factor of 2 smaller than those computed from the phase-diversity 
reconstructed images. 
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Fig. 10. Histogram of the CP at the corresponding positions of 
the Sunrise Ca ii H BPs in the photosphere. The dashed line 
marks the Scr noise level. The red, dot-dashed line corresponds 
to an exponential fit and the vertical solid-line indicates a mean 
values of 0.32% CP. 
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Fig. 11. Histogram of the unsigned longitudinal component of 
magnetic field, \Bco^iy)\, from the SIR inversion. The red, dot- 
dashed line represents an exponential fit to the right w^ing of the 
distribution, i.e., for B > 100 G. The vertical line indicates the 
mean field strength value. 



4.3.6. Longitudinal Magnetic Field 

As a final result, we perform an analysis of the Stokes pro- 
files using the SIR code to compute the unsigned longitudi- 
nal component of magnetic field, \Bco^{y)\ (see Sect. 13.41) . We 
show \Bco^{y)\ since this value is retrieved more reliably than 
B, whose values can be adversely aff'ected by noise in Stokes 
Q and U. However, we should keep in mind that the \Bco^{y)\ 
values may not be completely reliable, since the BPs either are 
not fully spatially resolved or the polarisation signals of the 
small BPs under study are too close to the noise level. The for- 
mer issue causes the inversion code to underestimate the true 
field strength. The second issue is caused by the sometimes ex- 
treme line weakening of the highly temperature and magneti- 
cally sensitive IMaX Fe I 5250.2 A line (Sheeley 1967; Stenflo 
1975; 'Shelv ag et al.l 120071: iLagg et al.l I20ia |Riethmuller_etalJ 
2010. 2012). The fact that such point-like magnetic features are 
bright, suggests that they are the cross-sections of kG magnetic 
elements whose field strengths are underestimated, very likely 
because their true diameters in the photosphere are below the 
Sunrise spatial resolution. A mor e thorough discussion is pro- 
vided bv lRiethmuUer et al.l (l2012l). 

Fig.[TT]depicts the distribution of \B cos(7)| for the BPs observed 
on 9 June 2009. This histogram shows a mean value of 142 G 
indicated by a vertical line in Fig. [TT] An exponential fit to the 
points with B > 100 G results in an e-folding width of 87 G. 
A similar conclusion as for the CP distribution is made 
w hen our longitudin al magnetic field is compared with that 
of Riethmiiller et al. 1 ('2012), who found stronger fields ranged 
between 25 - 1750 G and with a mean value of 544 G. 

4.4. A Relationship Between Intensity and Horizontal Velocity 

The only statistically significant correlation obtained from data 
of all BPs was found between the maximum values of intensity 
and maximum proper motion velocity of all identified BPs. 
The maximum intensity values of all identified BPs are plotted 
versus their maximum horizontal velocities in Fig. [121 The linear 
regression fit to the data (red solid line) shows an inverse corre- 
lation between maximum intensity and maximum proper motion 



velocity of the BPs with a correlation coeflScient of -0.49. The 
blue dot-dashed lines in this plot indicate the confidence bands 
with 95% confidence level which show an interval estimate for 
the entire regression line. Clearly, the faster a BP moves, the less 
bright it is. In particular, the combination of very bright and very 
mobile BPs is not present. 
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Fig. 12. Relationship between the maximum horizontal velocity 
and the maximum intensity of the identified BPs. The solid (red) 
line illustrates a linear regression fit to the data with 95% confi- 
dence bands (blue dot-dashed lines) around it. 



5. Kink Wave Excitation Due to High-velocity 
Excursions 

We assume that the migration of our Ca ii H BPs is a proxy for 
the motion of the cross-section of the underlying magnetic el- 
ements or flux tubes in the layer sampled by the Sunrise/SuFI 
Ca II H passband, i.e., the cross-sections of magnetic elements 
at around the height of the temperature minimum/low chromo- 
sphere. 

Furthermore, we found that all studied BPs have velocity vari- 
ations similar to the one plotted in Fig. O?. In particular, the 
maximum (peak) velocity observed for each BP is often much 
larger than its mean value and is reached for only a brief time. 
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This corresponds to the occurrence of intermittent pulses in the 
horizontal velocity. In Fig. [13] the mean and the peak values of 
horizontal velocity of all identified BPs are plotted versus their 
observed lifetimes. In both plots, higher velocities are mostly ob- 
served for shorter lived BPs. The peak horizontal velocity ranges 
between 2.6 and 15.5 km with a mean value of 6.2 km 
for all BPs. In addition, Fig. \T3h shows that almost all BPs 
move fast at least once during their lifetimes. The horizontal 
solid (red) and dot-dashed (blue) lines in both panels in Fig.[T3l 
indicate the rough limits distinguishing fast, intermediate and 
slow BP motions, with fast motions implying mean proper mo- 
tion speed ^ 3 km and intermediate motions between 2 and 
3 km . T hese limit s were chosen following .Choudhuri et al.l 
(Il993ah and lKalkofenl(ll997h . 
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Fig. 13. (a) The mean horizontal velocity, and (b) maximum 
(peak) horizontal velocity versus the observed lifetimes of the 
studied BPs. The horizontal solid (red) and dot-dashed (blue) 
lines indicate the approximate limits for fast and intermediate 
BP velocities, respectively (see main text). 



Note that vibration-induced image jitter introduced by instru- 
mental eff'ects (i.e., due to brief loss of pointing-lock of the tele- 
scope) into the Sunrise observations is much smaller than our 
observed jerky motions (Berkefeld et al. 20lT|). 
Rapid horizontal motions are of interest in particular if they oc- 
cur as a pulse, i.e., in combination with a rapid acceleration, 
since such motions have been proposed t o exc ite magnetoacous- 
tic kin k waves along flux tubes ( Spruitl ll98ll: iChoudhuri et al.l 
fl993ih . Such transversal kink waves transport energy and can 
trave l into the up per chromosphere without significant reflec- 
tion ("S pruitlll^SlI ). Although the reflection at the transition layer 
decreases the energy flux transported by the kink wave, it may 
be interesting for quiet coronal heating (Cho udhuri et al.l 1993b) . 
These early results have been given a strong boost by results of 
iDe Pontieu et al.l (l2QQ7l) that transverse magneto-hydrodynamic 
(IVLHD) waves along spicules, which may propagate along ex- 
tensions of BPs to higher layers, carry sufl S cient e nergy to heat 
the corona. According to IChoudhuri et al.l (Il993ab . proper mo- 
tions of BPs faster than ^ 2-3 km (happening within a time 



shorter than 100 sec) are particularly efl'ective in exciting kink 
waves. 

This mechanism can be imagined as the rapid jostling of flux 
tube by granules in the photosphere whic h impulsively excites 
oscillations in the tube ('H asan etaP 1200 0). Then, the oscilla- 
tions (waves) propagate upwards while their velocity ampli- 
tudes increase exponentially until they reach values compara- 
ble to the tube speed for kink waves in chromosp heric lay- 
ers, where the kink waves become nonlinear (e.g., iKalkofenI 
1997). In the middle of the chromosphere, these transver- 
sal kink w aves may couple to the l ongitudinal magnetoacous- 
tic waves (Ulmschneider et al. 1991) and dissipate by forming 
shocks (e.g., Zhugzhda et a l. 1995). 

IChoudhuri et al.l (Il993albl) carried out a theoretical analysis of 
the transversal kink waves induced by the flux tubes footpoints. 
They derived the following expression that gives a rough esti- 
mate of total energy flux (Fe) carried by the kink waves: 



6 5 X 10^6 ( PQ \ ( ^0 ] ( ^0 \' 
\ 10-7 g cm-^ I \ 105 km^)\l km s'^ I 



10 ^ gem ^ 



em ^ s ^ 



(6) 



where po is the atmospheric density, Aq is the cross sectional 
area of the flux tube (i.e., area of the BP), vq is the maximum 
velocity of a horizontal pulse, and H is the scale height. The pa- 
rameters n and / in Eq. (O represent the number density of BPs 
and the frequency of motions (i.e., pulses), respectively. /I is a 
dimensionless parameter related to the rapidity of a BP, i.e., a 
measure of the strength of the BP motion. This parameter is de- 
termined by /I = vo/(oJcL), where ojc is the cut ofl' frequency at 
a certain layer in the atmosphere and L is the BP displacement 
within the pulse duration. The dimensionless function F(A) is the 
asympt otic energy which depend s on the employed atmospheric 
model. IChoudhuri et al.l (Il993bl) examined three difl'erent mod- 
els and finally concluded that an isothermal atmosphere is a rea- 
sonable model for kink waves propagation in the solar atmo- 
sphere induced by rapid flux tube footpoint motions. Therefore, 
we employ the isothermal atmosphere model in our calculations. 
Since we are aiming for a rough estimate of the energy flux gen- 
erated by our high-velocity Ca ii H BPs, we are satisfied with 
approxim ate values of the various parameters entering Eq. 
Following I Wellstein et al.l (Ll998), we take the scale height H in 
an isothermal atmosphere of 7000 K (as an upper limit), equal 
to 150 km. At the height of formation of the Sunrise/SuFI 
Ca II H BPs (i.e., a height corresponding roughly to the tempera- 
ture mi nimum), we use the d ensity from the VAL-C atmospheric 
model (IVernazza et al.ll 19811) : po ~ 6 x 10"^ g em~^. 
For the acoustic cut off* frequency at the temperature minimum 
we use 3 x 10"^ ^"^ (Spruit 1981). 

Then, we analyse the horizontal velocity variations of each BP 
and look for significant pulses, i.e., when the proper motion ve- 
locity exceeds a certain limit for a maximum duration of 100 see. 
We compute the pulses occurring in 1 km s~^ wide velocity bins 
between values of 2 and 16 km s~^. The pulse parameters, i.e., 
the peak horizontal velocity and the area of each BP at that time, 
as well as the displacement which the BP suff'ers during each 
pulse, are recorded. The frequency of pulses for each bin is then 
computed as the ratio of the number of pulses in the bin and the 
integrated lifetime over all BPs. As mentioned earlier, the iden- 
tified BPs in this study have a number density of 0.03 (Mm)~^. 
The parameters vq, Aq and L are computed as the average over all 
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Table 4. Characteristics* of Ca ii H BPs pulses and estimated values of energy flux for difl'erent velocity bins. 



Velocity Bin 


Np 


vo 


h 


Ao 


L 


/ 


A 


F{A) 


Fe 


[km s~^] 




[km s~^] 


[s] 


[km^] 


[km] 


[mHz] 






[erg cm~^s~^] 


2-3 


252 


2.5 


42 


16000 


55 


4.00 


1.5 


0.33 


4.1 X 10^ 


3-4 


206 


3.5 


43 


16500 


61 


3.00 


1.9 


0.48 


5.9 X 10^ 


4-5 


135 


4.4 


43 


17100 


70 


2.00 


2.1 


0.54 


6.1 X 10^ 


5-6 


84 


5.4 


47 


17600 


70 


1.00 


2.6 


0.70 


5.0 X 10^ 


6-7 


31 


6.4 


42 


16900 


81 


0.50 


2.6 


0.71 


2.4 X 10^ 


7-8 


15 


7.5 


33 


18300 


72 


0.20 


3.4 


0.98 


1.4 X 10"^ 


8-9 


20 


8.4 


38 


18600 


80 


0.30 


3.5 


1.00 


2.4 X 10"^ 


9-10 


11 


9.5 


40 


18600 


74 


0.20 


4.3 


1.25 


1.4 X 10^ 


10-11 


7 


10.4 


22 


18600 


74 


0.10 


4.7 


1.39 


1.0 X 10^ 


11 - 12 


2 


11.4 


47 


25000 


46 


0.03 


8.2 


2.58 


2.7 X 10^ 


12-13 


2 


12.4 


28 


19800 


70 


0.03 


5.9 


1.81 


3.5 X 10^ 


14-15 


3 


14.2 


49 


19700 


76 


0.04 


6.2 


1.90 


6.5 X 10^ 


15-16 


2 


15.7 


33 


22300 


65 


0.03 


8.0 


2.49 


4.8 X 10^ 


2-16 


770 


4.1 


42 


16800 


63 








3.1 X 10^ 



* Np'. Number of pulses; vq: Mean horizontal velocity; tp'. Mean lifetime; Aq: Mean area; L: Mean displacement; /: Frequency of pulse occur- 
rence; A: Rapidity parameter; F(A): Asymptotic energy; Fe'. Net flux 



the peak horizontal velocities, all the areas and all the displace- 
ments determined for all pulses in all BPs within each velocity 
bin, respectively. We extrapolate the plot of F(A) versus A from 
Fig. 6 in iChoudhuri et al. (J[993a) to obtain the value of F(A) for 
the computed As. The values of the obtained parameters for each 
velocity bin are presented in Tabled The integrated energy flux 
is then computed over all velocity bins larger than a given speed. 
We do the analysis for two speed limits: one for v = 3 km 
and one for v = 2 km representing fast and intermediate B P 
motions, respectively (IChoudhuri et al.|[T993al:lKalkofenlll997l) . 
We find an integrated value of Fe ~ 2.7 x 10^ erg cm~^s~^ (or 
^ 270 W m~^) for pulses with v ^ 3 km and a marginally 
larger value of 3.1 x 10^ erg cm~^s~^ {ox ^ 310 W m"^) 
for pulses with v ^ 2 km s~^. These estimated values of total 
energy flux are compatible with the estimated energy flux den- 
sity of 1 - 3 X 10^ ergs cm'^s'^ (100 - 300 W m"^) required to 
heat the quiet corona an d/or drive the solar wind ( Hollweg 1990; 
iHansteen & Leej|l995h . Our value is a lower limit since we only 
considered small, internetwork magnetic elements. 
Table |4] is in agreement with the finding of Choudhuri et al. 
(1993b), that faster motions, though more infrequent, provide 
more substantial contributions compared to the slower ones (cf. 
Eq.O. 

Note that if the formation height of the Ca ii H channel were 
to lie lower, e.g., in the mid photosphere where the atmospheric 
density is much higher, po ~ 5x10"^ g cm~^, then the energy flux 
estimated based on Eq. © would be about one order of magni- 
tude larger. 

iMuller et"aD (1 19941) reported a very much larger value of Fe = 
2000 W m~^ for their ph oto spheric network BPs with v > 
2 km whereas IWellstein et al.l (11998 ) found a value of Fe = 
440 W m~^ for their K2y BPs with mean horizontal velocity of 
6 km s~^. The latter value is highly afl'ected by the formation 
height assumption of 1100 km for their Ca ii K passband, and 
may not be relevant in any case since it is unclear w hich fractio n 
of their features really is magnetic. The results of iMuUer et al.l 
(.1994) suggest that it would be worthwhile to follow also larger 
BPs as found in the network in Sunrise data. 



6. Conclusions 

In this paper, we have provided observational properties of in- 
trinsically magnetic and highly dynamic small bright point-like 
features in the Sunrise Ca ii H passband. BPs of the type studied 
here, i.e., with diameters below 0.3 arc sec, are present mainly in 
the internetwork. We have, for simplicity, further restricted our 
analysis to those BPs that do not merge or split in the course of 
their observed lifetime. They cover only 5% of the Sun's surface 
and radiate on average 1.48 times more than the average quiet- 
Sun in the 1.8 A broad Sunrise Ca ii H filter, so that 7.5% of the 
radiative energy losses in lower chromospheric layers are due 
to such features (assuming Ca ii H to be representative of other 
spectral lines formed in that layer). This is negligible. However, 
motions of these features may be responsible for kink waves 
travelling into the corona, and for the braiding of field lines lead- 
ing to microflares (Parker 1988), so that studying them can po- 
tentially help us in our understanding of chromospheric and pos- 
sibly coronal heating mechanisms. 

We applied a set of stringent criteria in order to separate the 
(magnetic) BPs we wish to study from all other bright, point- 
like features. In particular, we only consider BPs smaller than 
0.3 arcsec. 

This limited us to a smaller number of BPs, but ensured that they 
are not oscillations or wave-like features, nor the product of re- 
versed granulation. 

Nearly all the BPs for which we could test this were found to 
be associated with signi ficant photo spheric ni agnetic field. This 
result agrees with that of Ide Wij n et al.l (l2QQ5l) , who showed, al- 
though at lower spatial resolution, that internetwork Ca ii H BPs 
were located inside magnetic patches. They also found a good 
spatial coincidence between their Ca ii H BP patches and the 
ones they fo und in G-band images. T his is also the case in the 
Sunrise data ('Riethmii ller et al .1120101) . 

These small magnetic BPs in internetwork areas have long life- 
times with a mean value of 673 sec, much longer than any val- 
ues in the literature for Ca ii H BPs, but in good agreement with 
the mean lifetime of about 9 min for photospheric internetwork 
BPs (Mullerl983) as weU as a mean lifetime of 10 min for the 
internetwork magnetic elements (IME; de Wijn et al. 2008), of 
which only 20% are associated with G-band and Ca ii H BPs. 
Our internetwork features appear to live less long than (photo- 
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spher ic) network BPs which last on average 1080 sec (iMuUed 
11983"), however, a comparable mean lifetime of 9.2 min for net- 
work GBPs was found by Nisenson et al. ( 2003). In addition, we 
showed that the Ca ii H BPs move horizontally with an average 
speed oil.lkm . We suspect that the muc h hig her horizontal 
velocities reported by IWellstein et al.l (Il998l) and ISteffens et al.l 
(1996) refer to a different type of feature, which is probably non- 
magnetic. Furthermore, an inverse correlation between the maxi- 
mum values of intensity and horizontal velocity of the identified 
BPs was observed, i.e., the BPs are brightest when they are at 
rest. 

It was shown that in addition to the considerable mean values of 
the horizontal velocity of the subsonic BPs, almost all Sunrise 
Ca II H BPs move with a much larger peak horizontal veloc- 
ity (up to 15.5 km s~^) for at least a short period of time. On 
average a BP spends 3.5% of its life travelling at a supersonic 
speed. These motions are not regular but often correspond to 
pulses of rapid horizontal motion. Such jerky, pulse-like motions 
eflSciently give rise to kink waves travelling al ong a magnetic 
flux tube according to IChoudhuri et al.l (Il993allb ). They calcu- 
lated that the efficiency of kink wave excitation depends on the 
rapidity of the BP motion with respect to the cutoff' frequency 
of the atmosphere. We computed a rough estimate of the en- 
ergy ffux generated by the rapid BP motions (i.e., for pulses with 
V > 2km s-^) of ^ 3.1 x 10^ erg cm-^s-^ (or ^ 310 W m'^), 
which is sufficient to heat the quiet corona (energy ffux density of 
1 - 3x 10^ er^s cm-^s-^ or 100-300 W m"^:IWithbroe & NovesI 
119771: iHollwegll 19901: iHansteen & Leedl 19951) if the kink pulses 
propagate into the corona and dissipate their energy there. Our 
estimated value is an order of magnitude lowe r than the energy 
flux tr ansported by Alfven waves estimated by lDe Pon tieu et al. 
(l2007l) from observations of type ii spicules (4 to 7 kW m~^). 
However, the mechanisms by whi ch these diff'erent energy 
fluxes are generated may be diff'erent. [Hasan & van BallegooijenI 
(l2QQ8h argued that the heating in Ca ii H BPs, caused by weak 
shocks occurring at short time intervals (less than 100 sec) in 
magnetic ffux tubes, is in contrast with the long-period waves 
which are considered as the spicules driver ( De Pontieu et al.l 
|2004|) . Also, by restricting ourselves to small BPs, we may be 
missing the magnetic features carrying the most energy into the 
upper atmosphere. 

The mean value of the unsigned longitudinal component of mag- 
netic ffeld, \Bcos(y)l in the BPs was found to be 142 G, reach- 
ing 500 G only rarely. These values may be inffuen ced by some 
of the BPs not being fully spatially resolved (iRiethmiiller et al.l 
l2012h and/or the polarisation signals lying close to the noise 
level at some times, as can be seen in Fig.fTOl 
All distributions of the parameters studied here, except that of 
diameter, display an exponential fall-off' towards larger values, 
implying random stochastic processes. The shape of the size his- 
togram is partly determined by our restriction to features smaller 
than 0.3 arc sec, and may well be aff'ected by the fact that many 
of these features could be unresolved in the photosphere, as is 
suggested by their weak ffelds (see Riethmiiller et al. 20 12). 
We identiffed 7 persistent ffashers (iBrandt et al.l Il992b in our 
sample of BPs. We ffnd that the ffashers are normal BPs that 
diff'er mainly in that they have smaller average values of inten- 
sity, so that their brightness drops below the threshold at some 
moments in time and they seem to disappear temporarily. 
We expect that granular and intergranular motions are primar- 
ily responsible for the horizontal motions of the corresponding 
ffux tubes rooted in intergranular lanes. Interestingly, the life- 
time distribution of our Ca ii H BPs (range between 16 7 and 
1507 sec) is comparable with that found by Hirzberger et al.l 



(Il999l) for granules, who showed a range of 168 - 1800 sec for 
granular lifetimes with only a relatively small number of long- 
lived granules with lifetimes longer than 1200 sec. However, 
they found a mean lifetime of about 360 sec for the granules 
which is smaller than the lifetime of our Ca ii H BPs. In addi- 
tion, the existence of horizontal convective supersonic ffows at 
the boundaries between granules and intergranular lanes ( e.g. , 
Cattaneo et a l. 198 9; Solanki et al. 1996; Nordlund et al. 20091: 
iBellot Rubioll2009l: Ivkas et al.ll2()riD may explain our observa- 
tions of fast motions BPs and supersonic pulses which can be 
due to the impact of fast granular ffows on the ffux tubes. 
Summarising, we studied the structure and dynamics of the 
smallest currently observable magnetic bright features at the 
height sampled by the Sunrise/SuFI Ca ii H passband (i.e., a 
height corresponding roughly to the temperature minimum). We 
determined and discussed the statistical distributions of the BPs' 
properties including intensity, horizontal velocity, size, lifetime, 
polarisation signals as well as magnetic ffeld strength. We found 
an anti-correlation between the maximum proper motion veloc- 
ity and the peak intensity values of the BPs. With the help of an 
advanced image processing technique, selecting the actual mag- 
netic BPs based on stringent criteria and using data unaff'ected 
by seeing from Sunrise, we provided accurate horizontal veloc- 
ity proffles of the Ca ii H BPs. These proffles revealed indications 
of fast pulses (i.e., rapid BP motions on short time scales) which 
may contribute to coronal heating by exciting kink waves in the 
corresponding ffux tube s. The e nergy ffux that we estimate fol- 
lowing .Choudhuri et a l] (ll993afb ) is marginally sufficient to heat 
the quiet corona. Note that we restricted ourselves to small inter- 
network magnetic elements, so that we are probably including 
only a small fraction of the energy ffux in our estimate. Hence, 
in addition to other possibilities of coron al heating, e .g., braid- 
ing of ffeld lines leading to microffares ( Parkedl 1988b . the jerky 
motions of the BPs studied in this work, can excite waves that 
potentially carry enough energy to contribute to the heating of 
the quiet corona. Future work will aim at detecting such waves. 
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